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We quantitatively analyze the meson mass inequality relations of two dimensional gauged four 
fermi models in the large A*' limit. The class of models we study includes the 't Hooft model, the 
chiral and non-chiral Gross-Neveu models as special points in the space of field theories. Cases 
where the chiral symmetry is spontaneously or explicitly broken are both studied. We study the 
meson mass inequality quantitatively and define a susceptibility which allows us to systematically 
analyze the inequality. In the generalized Gross-Neveu model limit, we derive an analytic expression 
for this susceptibility. Even though no analytic proof of the validity of the classic mass inequality 
exists for the generic case, the mass inequality is found to be positive throughout most of the 
^~>) ■ parameter space. We point out that the inequality might be negative in certain cases. 

(N : 
< ■ 

Determining the properties of composite particles in relativistic field theories, such as mesons or baryons in QCD, 
is an essentially non-perturbative problem and is quite difficult to address from first principles. While impressive 
■ progress has been made in numerical approaches to the problem it is highly desirable to also have a more 
I , analytic understanding of the properties of the composite particles. An important relation that can be applied to 
<^ composite particles in vector like theories, like massive QCD, are the mass inequality relations While these elegant 
, inequalities are quite useful and have therefore been well studied (for a recent review, see, |^), there seems to be little 
' understanding regarding their quantitative behavior. Furthermore, while the proof of inequalities does not apply to 
" chiral, explicitly left-right asymmetric theories or theories with Yukawa couplings, no examples of theories wherein 
. the inequality has been shown to be negative exists amongst relativistic quantum field theories. The mass inequalities 
^-^ ' provide non-trivial insight into the dynamics of the interacting model. An understanding of the behavior of the 
inequalities will undoubtedly further our understanding of the spectrum of bound states in relativistic quantum field 
theories. 

Whether the inequalities might or might not be positive in chiral models is of import, in particular to supersymmetric 
theories. In supersymmetric theories, one often 'solves' the gauge hierarchy problem by invoking the chiral symmetry 
O il of the fermions, which is related in turn to the mass of the scalar bosons by supersymmetry. While some non- 
D ' perturbative aspects of supersymmetric gauge theories are recently being clarified relations analogous to the 

mass inequality relations seems not to be known. Such relations, if they exist, should shed light on the properties 
^ , of the spectrum on supersymmetric theories and on the possibility of spontaneous breaking of supersymmetry. In 
attempting to extend the mass inequalities to supersymmetric theories, we believe a deeper understanding of its 
. , properties to be crucial. Of course, non-supersymmetric theories are of interest on its own, one of the reasons being 
> that the low energy world is not supersymmetric. 

While, needless to say, computing the mass inequalities in full QCD would be of import, not to mention very 
interesting, this is a daunting task. One approach is to compute the mass inequalities and develop an understanding 
of their behavior in analytically solvable relativistic quantum field theory models, such as the classic 't Hooft model 
1^ and the Gross-Neveu models [Q. This is the approach we shall adopt in this work. These models are tractable 
yet non-trivial and have proven to be quite instructive by providing physics insight into the non-perturbative aspects 
of field theories . Dynamics of gauge theories is certainly of import especially since it is an integral part of the 
Standard Model. ^Iso, theories of four fermi interactions have been serving an important role in particle physics and 
other fields of physics |jll|,|l^ . 

In this work, we shall analyze the properties of mass inequalities of gauged four fermi models in (1 -t- l)-dimensions, 
using the large N limit. In these models, the properties of the "meson" states can be reduced analytically to the 
problem of solving mathematical equations ||l^,^ . This will allow us to study the problem analytically, even if the 
final equations need to be solved numerically. The parameters in this family of models which we can arbitrary control 
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are the scalar and the pseudo-scalar four fermi couplings, the gauge coupling and the "quark" masses. The class of 
models we study contains the 't Hooft model ||] and the chiral and non-chiral Gross-Neveu models for particular 
choices of the parameters. It includes cases where the chiral symmetry is spontaneously or explicitly broken. This is 
of particular interest, since the Gross-Neveu model is known to be equivalent to a Yukawa model where the scalar 
field develops an expectation value ||^. This is precisely the kind of situation where one might suspect that the mass 
inequality may be violated |^,^. 

In order to study the mass inequalities quantitatively and systematically, we need a measure for the size of the 
inequality. We shall find such a measure in a parameter which we shall call the 'meson mass susceptibility' to be 
explained below. This parameter will allow us to compare the inequalities of various field theories in a natural manner. 

The work is organized as follows: In section ||, we shall define and explain the meson mass susceptibility. In 
section III, we discuss mass inequalities in quantum mechanics. In section [V, we first treat the simpler case of 



generalized Gross-Neveu models and analyze the mass inequalities there. In particular, an analytic form of the mass 
susceptibility will be explicitly presented. In section we work with the general gauged four fermi models. We 
first summarize the class of models we study and explain how the meson spectrum can be computed. Particular care 
has been paid to presenting the methods we use explicitly for further possible use. We then compute and study the 
properties of the mass differences. We end with discussions of the results in section A short appendix on some of 
the technical aspects of the computation is included. 



II. THE MESON MASS SUSCEPTIBILITY 



In studying the mass inequalities quantitatively, we need a quantitative measure of how "large" the inequality is, 
in order to compare within the field theory space. This encodes information, intuitively speaking, on how strong the 
attractive interactions in the theory are. The mass inequality we consider is 

r _ ifJ-aa + ^J'bb) n 
0^^ab = Mab 2 

which is known to be positive for vector-like gauge theories Here, denoting the constituents as q's, fiab is the mass 
of the lightest meson that overlaps with the q^qt state, and so on. 

This quantity is dimensionful and depends not only on the model, but also on the difference of the masses of the 
constituents. The quantity may be made dimensionless trivially by taking the ratio of the inequality with a meson 
mass. However, this inequality may become large just because the constituent mass difference is large, so is not the 
most appropriate parameter for investigating the intrinsic dynamics of the theory. In fact, if we consider the field 
theory space to be parameterized by the couplings of the model which include the masses, the mass inequality is not 
a local quantity in the parameter space. It is more natural to define a local parameter in the field theory space. Let 
us define the mass squared of the constituents to be 

= 1/2(1 -I- A), = M^{1 - A) (2) 

The meson mass difference Sfiab is even under the interchange of Ma and M}, and is therefore an even function of A. 

We shall characterize the inequality by a parameter we refer to as the "meson mass susceptibility" . This quantity 
is defined by, 

,. ^fJ-ab '^fJ'ab - ifJ-aa + fJ-bb) 

TZ = hm —5 = hm —r- , (3) 

which is a function of the couplings and the mass of the constituent quark. The susceptibility we defined above is 
also useful from a practical point of view: Since the mass inequality is an expansion in the mass difference squared, 
the susceptibility together with the meson mass for equal mass quarks, distill the meson mass information when the 
quark mass difference is not too large. 

We should perhaps here discuss the relation of this susceptibility to the global properties of the mass inequality, 
namely when the mass differences are arbitrary. A natural question is whether the positivity of the susceptibility 
in a parameter region guarantees the positivity of the inequality when the mass differences are large. In quantum 
mechanics, the situation is quite simple; if the susceptibility is positive everywhere, the mass inequality is positive for 
arbitrary mass differences. This can be derived from the convexity of the meson mass with respect to the reduced 
mass of the two quarks. In quantum field theory, however, no such argument exists in general, since the meson mass 
needs not and will not depend only on the reduced mass of the two quarks. To make an analogous argument in 
quantum field theory, we need further information regarding the relation between the meson mass and the quark 
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masses. While it seems quite natural to assume that relations exist such that the positivity of the local susceptibility 
guarantees the positivity of the mass inequality globally, we do not know if this is in fact true. In practice, we have 
found no counterexamples to this statement. 



III. MASS INEQUALITIES IN QUANTUM MECHANICS 

In this section, we briefly discuss mass inequalities in quantum mechanics (see also |^). While the discussion is not 
necessary for computing mass inequalities in relativistic field theories, we feel that it is nonetheless quite instructive 
and provides a broader perspective on mass inequalities in quantum theories. Also, the mass inequalities in relativistic 
field theories should reduce to that of quantum mechanics in the non-relativistic limit. As such, some of the results 
here will be later compared to those from the full quantum field theory below. It should be noted, however, that 
phenomena such as symmetry breaking which plays a large role in (1 + l)-dimensional gauge theories studied in this 
work are essentially quantum field theoretical so that quantum mechanical behavior is not sufficient for understanding 
the full relativistic behavior, even qualitatively. 

In quantum mechanics, the problem of two body bound states under a local potential reduces to a model with the 
hamiltonian 

where M12 denotes the reduced mass, l/A/12 = 1/Afi + \/M2- We will analyze one dimensional models, but similar 
analysis can be applied to higher dimensional models. 



1. Infinitely deep square well potential 

The potential of the model is 

The spectrum of the bound states is known to be i?i2,n — 2M ^ {n = 1,2, . . .). This is somewhat trivial but an 
interesting case. The susceptibility TZ — and this we can uncferstand as the signature of the model being free within 
the well. 



2. Delta function potential 

The delta function potential 

V{x) = -Vo5{x) {Vo > 0) (6) 
has a bound state with the binding energy — Mi2Vo^/(2?l^). 

. ^ Eaa+ Ebb [Ma - Mbf 

d^lab = -C^ah -Z = —79-^77 , . ,r \ \' I 

2 {Ma + Mb) 



This leads to the susceptibility 

32(1-^078) 



The susceptibility increases with larger Vq, as expected. In the non-relativistic limit, Vq 1. 
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3. Monomial potentials 



Let us also discuss potentials whose behavior is governed by a monomial 

V{x)=Ax^, A7>0 (9) 

7 needs not be an integer but 7 > —2 needs to be satisfied for sensible physics behavior. A-f > needs to be imposed 
for the existence of bound states. 7 = 2 and —1 corresponds to the harmonic oscillator and the three dimensional 
Coulomb case, respectively. 

We can use the uncertainty principle to crudely estimate the bound state energy as 



We can obtain the susceptibility from this energy as 



h,-, + 2 /jA\-> + '^ 2(^ + 1) 

7^~— — p^rr Kr M"^+^ > o (ii) 



8(7 + 2)c2 V 2 

in the non-relativistic limit. While the derivation is not rigorous, in the harmonic oscillator and the Coulomb cases, 
the susceptibilities agrees with those obtained from exact methods. 



IV. GENERALIZED GROSS-NEVEU MODELS 



In this section, we analyze the mass inequalities in the generalized Gross-Neveu models, described by the lagrangian 

Nf 2 2 

^ = ^^}ii^-^f)^f + ^ XI (V'/'V'/)(V'/V'/') - y 51 (V'/'75V'/)(V'/75V'/') (12) 

/=1 /,/' = ! = l 

In addition to the flavor indices /, /' denoted explicitly in the above formula, the fermions carry an additional internal 
space index, the 'color' index (l,2,...iV) which has been suppressed in the notation. This index should not be 
confused with the flavor index. We take the large N limit while keeping a^N^a^N fixed. When mf — Q,a\ — 
the model reduces to the original Gross-Neveu model and when to/ = 0,a^ = a|, the model reduces to the chiral 
Gross-Neveu model with continuous chiral symmetry. When to/ = and the couplings are not equal, we are left with 
discrete chiral symmetry in the model. We need to consider multiple flavors for the analysis of the mass inequalities. 

This class of models is included in the gauged four fermi models we deal with below and the analytic methods 
discussed there can be applied here also. However, the generalized Gross-Neveu models can be solved completely 
analytically using different methods than the gauged four fermi model case, so we shall discuss it separately. Here, 
we shall need the spectrum in the general case when two flavors have different masses, m\ ^ m\, and a? ^05, which 
was not solved explicitly in |l^]. We shall present the spectrum and analyze the mass inequalities. 

Let us consider a meson bound state of constituents with masses, Mi,M2. These constituent masses are physical 
fermion masses that include the effects of spontaneous chiral symmetry breaking that occurs dynamically in the 
Gross-Neveu model. We dispense with the derivation here, but the Bethe-Salpeter equation for the meson state can 
be solved algebraically to obtain the meson "wave function" , (p{x) as 

(p(a;) = v'f") +<^(i)(l-2a;) + (^(a;), (0 < a: < 1) (13) 
_ + y(i)(l - 2x)) + 2 {Ml - M|) yd) 

^i-Z X 1 — X 

where </)(°\ (/j^-*^) are constants and (p{x)/ [x{l — x)] is integrable at a; = 0, 1. The meson wave function satisfies the 
following boundary conditions 

b_ (1 + 4G)6+ j 1,(^(1) j -y/^x(l-:r) G)[b^ b+)[l-2x) ^^^) 
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Here, we used the notation G = a?N/{AT:), G5 = a|A'"/(47r) for the renormalized couplings and defined b± 



(Ml ± M2) /2. When the couphng constants are equal, a 



^5; 



or when the masses are equal. Mi = M2, the 



two boundary condition equations simply decouple, but do not in the general case. 
The boundary conditions lead to a secular equation 



det 



J12 — 2 



1 

G5 



1 

G5 



4Mi Mo 



■J12 +^T^m^+[-g^-g) 



J12 + 



In 



4M, Mo 



4M1M2 



Here, we defined 



J12 



/i^ dx 



-l-j?x{l 



M2(l 



+ M^x 



(16) 



(17) 



It should be noted that since the couplings G, G5 are dimensionless, the overall mass scale M can always be scaled 
out of the problem and only the relative masses have a physical meaning. The physical parameters of this quantum 
field theory are the two dimensionless renormalized couplings G, G5 and the mass ratio A/1/M2. 

Before we analyze the behavior of the mass inequalities, we first need to understand the behavior of the spectrum 
when the masses of the constituents are the same. In this case, the secular equation (|l6|) splits into two independent 
equations for the pseudo-scalar and scalar bound states, x ^-nd cr: 



cr : 



^4(M/mJ2 



■■ tan 



1 _ rh^ (i^JMf 

G ~ Jo l-fn^/MV-i 



(18) 
(19) 



It should be noted that only G5 (G) appears in the equation for x (f)- 

a exists as a non-tachyonic bound state only for G < —1/4. It is not clear whether the theory is unitary for negative 
G and we shall consider the region G > 0, so we shall not have much more to say on a. The original Gross-Neveu 
model corresponds to G — > —00 in our scheme and in this limit, /i^ — > 4M2. 

X exists as a bound state for any G5 > and < /z^ < 4M2. The dependence of the bound state mass on 
the coupling is plotted in Fig. |l|. This is the only bound state in the model for G, G5 > and corresponds to the 
Nambu-Goldstone like particle when the constituent masses are zero |0,|l^, as in the chiral Gross-Neveu model. It 
is the dependence of this meson state on the constituent masses that we shall investigate. As a side note, in a region 
we shall not investigate, there is an intriguing possibility when G5 > and G < —1/4, in some cases, the x uiass 
can be larger than the a mass. We do not know whether this can be achieved in a physically consistent situation. 
Another comment is perhaps appropriate; in the literature, the Gross-Neveu model (G ^ 0, G5 = 0) is often used as 
a prototypical simple model with a bound state. However, the original Gross-Neveu model has no binding energy for 
the meson and has barely a bound state. It seems to us that in fact, the simplest theory that may be considered in 
this family that is useful in analyzing bound state dynamics is G5 7^ 0, G = case. In this case, we have a bound 
state whose mass depends on the coupling as in Fig. |l|. 




FIG. 1. The behavior of the x meson mass, p-^/M with respect to the coupling, Gs. 
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The meson mass susceptibility may be obtained by perturbing the equation ( [T^ ) in the mass difference parameter 
A in Eq. (m). After some computation, we derive 



n 




i'^ i) Gk + 4G 



(20) 



We defined ( = Al'^/fiQ, where fig is the mass of the meson in the unperturbed case, when A = 0. The susceptibihty 
is independent of the overall mass scale M, since it can be scaled out of the problem. The susceptibility may be 
shown analytically to be positive for any G, G5 > 0. Since G, G5 are scalar and pseudoscalar couplings that can take 
on arbitrary values, the standard proof of the mass inequality does not apply to the models we study, except at 
special points. We believe that an analytic expression has not been derived for the mass inequality previously in any 
relativistic quantum field theory. 

It is interesting to check the asymptotic behavior of the susceptibility for small and large couplings. For small G5 
couplings. 



n ^ f3. [1 _ 8G5(1 + 4G) + 0{Gl)\ 



(21) 



This behavior is consistent with that for the 5 function problem discussed in Sect. 2. For large G5 couplings. 



7^ = 



Gs 



1 



2(4G+1) 24(4G+1)2 



o g; 



(22) 



The behavior of the susceptibihty with respect to G5 is shown for G — 0,0.1, 1, 10 in Fig. |[ The dependence on 
G is not strong; this is because the properties of the bound state x is governed mostly by the pseudo-scalar coupling 
G5. The crossover from G\ behavior to G5 behavior in the susceptibility can be clearly seen in the plot. 



100 



1e-06 



TT — I 1 — n — I 1 — n — I 1 — n — i 1 — r 




0.001 



100 



FIG. 2. The behavior of the meson mass susceptibility TZ with respect to the coupling G5 for the generalized Gross-Neveu 
models. The lines represent, from top to bottom, TZ for G = 0,0.1, 1, 10 respectively. 



V. GAUGED FOUR FERMI MODELS 



A. The model 



Let us now discuss the most general gauged four fermi model described by the lagrangian 

£ = -- tr (F^^F^^n +J2^f{il/)-mf)^f + y (i^^f)i^^f') " y E i^^5^f)i^l5^f') (23) 
/=i fj'=i f.f'=i 

We have gauged the internal index in the generalized Gross-Neveu model (^ so that when we set the gauge coupling 
to zero, we recover the generalized Gross-Neveu model discussed in the previous section. When we set G = G5 = 0, 
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we recover the 't Hooft model. We take the large- A'' limit in a manner similar to that of the previous section but also 
for the gauge coupling; namely, we keep g^N, a?N, a\N fixed while we take N to infinity. 

We again split the meson wave function as in ( |T3| ) . Then the Bethe-Salpeter equations for the meson bound states 
with the fermion constituents with masses Mi_2 may be obtained in a simple closed form |16|], 



^'^ip{x) = Hip{x) 



^ 132-1 

X 1 — X 



2ip 



(1) 



+ /32 + In 



l-x 



(24) 



satisfying the boundary conditions (p^. Since the gauge coupling g has the dimensions of the mass in (1 + 1)- 
dimensions, we introduced the dimensionless mass parameters j3i 2 = '"'-^i 2/(.9^-^)- "^^^ renormalized couplings G, Gs 
are defined in the same way as in the previous section. We should point out that all the parameters in this equation 
are finite renormalized parameters, so that the problem has been reduced to solving a somewhat complicated integral 
equation. The physical parameters in this theory are the three dimensionless renormalized parameters /3, G, G5 > 0. 

For later purposes, we also derive the matrix elements of the "hamiltonian" , iJ , in the most general case, when the 
couplings and the masses are arbitrary: 



1 



1 



4 VG ■ GsV ' 2 
G G5 



v,(o) v(°) + + (^(«) v^^) + ^« V(°)) 



dx 



Pi 



1 



/32 



l-x 



^+2^ g-g; 



dx dy 



(p'{x)ip{x) — P 



qJo [x-yf 



Lp'{x)(p{y) 



+ I dx2[ +/32 + ln 
'0 \ X 



(25) 



B. Methods for obtaining the spectrum 



In generalized Gross-Neveu models, the spectrum could be obtained by just solving an ordinary equation, albeit 
an transcendental one. In contrast, for the gauged four fermi model, we need to solve an integral equation which is 
technically more involved. Of course, this is to be expected, since the usual 't Hooft model, which is a simpler model. 



is solved in terms of an integral equation. To solve the integral equation (|24|), we employ two methods to be explained 
in this subsection, generalizing the methods used previously in the 't Hooft model p8| , p^ , p^ 16 1. With either of the two 
methods, we can solve for the spectrum and the wavefunctions of any of the meson states for arbitrary combinations 
of masses and couplings in the gauged four fermi models. By using the two different methods simultaneously, we are 
able to obtain a better control over the error in the results which inevitably arise when we solve the integral equation 
numerically, in addition to checking the internal consistency. We will be succinct and summarize the results. Even 
though the basic ideas are the same as those in the results are substantially more complicated since we need to 
treat the most general case which was previously not necessary. 

1. Variational method 

One method for solving the Bethe-Salpeter Eq. (^^, familiar from solving the Schrodinger equation, is the varia- 
tional method. We choose the basis functions {(pj\ j — 2,3, . . .} as 



'P2k{x) Cii + C2l(l - 2x) + 



B{k,k) 



. „^ (2fc + !)(!- 2a;) [a;(l-j:)l'' , 
(^2fc+i(a^) =ci2+C22(l-2x) + ^ ^1 ^ ^ (fc = l,2...) (26) 

Cy 's need to be determined to satisfy the boundary conditions ( |l5| ) as 
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Cll Ci2 
C2I C22 



5_ (1 + 4G)6+ Gy) 6. 

1 / -\{G - Gr^m + (32) + + iGsG) 

d\ i(G-G5)(/3i-/32) i(G-G5)(/3i+/32) + i(G5 + G)V^ 



(27) 



^l{G~G,m-p2) 



where 



(28) 



d = (1 + AG)bl - (1 + 4G5)6^ = [(G - G5)(/?i + /?2) + (1 + 2(G5 + G)) 
In the variational method, the problem of obtaining the meson states is reduced to solving an eigenvalue problem: 
ifi^Nki-Hki)wi=0, Hki = iipk,Hipi), Nki = {vk,w) kj = 2,3,4... (29) 

We will approximate the solution by using basis elements up to a certain number and check the convergence by varying 
the dimension of this basis space. With some work, the matrix elements can be computed to be 



I 



k + l B{k + l,k + l) 



^'''^^ " 3 2 V2fc + 1 ' 21 + 1 J ' 2{2k + 2l + l) B{k,k)B{l,l) 

_ 2 . 42 . C22 [ k I \ {k + l){2k + l)(2l + l) B{k + l,k + l) 

2k+iai+i ci2 3 2 ^2fc + 3 21 + i) 2{2k + 2l + l){2k + 2l + 3) B{k,k)B{lJ) 



1 k I 

N2k,2i+i = N2i+,,2k = cnci2 + 3C21C22 + 2(2fc + l)"i2 + 2(2[T3) 



C21 



(30) 



2k,2l 



-,VW2 



-(G + G5 + 8GG5)(/3i + /32) - -(G - G5) 



2c^2 



+ (/3i - /32)ci2 2cn + 



2k + 1 21 + 1 



(3i + p2 _ ^ B{k + l,k + l) 



kl 



B{k,k)B{l,l) 2{k + l) 



H2k+i,2i+i — -jy+hfh 



+ 2C22 [C22 - iPl - ^2)012] 



i(G + G5)(/3i + (32) + i(G - G5) 

k I \ f(3i+(32 \{2k + l){2l + l)B{k + l,k + l) kl{2k+l){2l + l) 

fcTl^r+T/^V 2 J 2k + 21 + 1 B{k,k)B{l,l) ^ 2{k + l){k + l + l) 



(31) 



-C22 



H2k,2l+1 - H2l + 1^2k - ^ - 02){G + G5) - 2C12 (C22 " (/?! " (32)Cl2) - ^12^ 

(R (2l + i)B{k + l,k + l) 



'2k + l ' 2'"'" {2k + 2l + l)B{k,k)B{l,l) 

When (3i = P2, the even and the odd sectors completely decouple. 



2. Multhopp's method 

Rather than using a variational method, we can expand the meson wavefunction and solve the eigenvalue problem 
directly [ p^ . Defining x = (l + cos^)/2, the wave function can be expanded in a manner consistent with the boundary 
conditions as 

CK K \ / K K \ K 

Cll ^ - C12 ^ u„ j - 27r I C21 ^ - C22 ^ w„ j COS 6* + ^ w„ sinn6' (32) 
n:odd n:cvcn / \ n:odd n:cven / n— 1 

where Cij's were defined in (p7|). This reduces the Bethe-Salpeter equation (^) to 

K 

J2 [f^^Pn{0) - M,^{d)] Vn=0 (33) 

where 
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P„(6') EE smn^ + 27r <^ ii ^34^ 
^ ' L — ci2 + C22 cos n: even 

~ / /?! — 1 /32 — 1 \ ^ n sin 710 „ , l + cos6'\ r — C91 

M,,ie) = 2[-tl - + -tl ^ sinn^? + 2^— — + 47r /3i-/32 + ln- x /^i 

V 1 + cos 1 — cos & J smO \ 1 — cos J L C22 



COS0 n: odd 
— C12 + C22COS0 n: even 
/3i - 1 , /32 - 1 \ . „ , ^_nsinn6' , f ^ ^ , i„ 1 + cos6'\ f -C21 n: odd 

n: even 

The above equation (^) is still a functional equation, with the dependence on the parameter 9. 
This can be further reduced to a generalized matrix eigenvalue problem 

{fi^P ~ M) v^O (35) 

The matrices are defined as 

K K 
Pmn=Y.9ra{ei)Pn{ei), Mr^n = 9m{0l) MM , 9,=Trj^ (36) 

1=1 1=1 

The function gm(0) is arbitrary, but using functions with the property gm{9i) = (— l)™^^.9m(0A'+i-;) simplifies the 
matrix elements. With this condition, the matrix elements are 

'' cii (to, n: odd) 

C22COS0; (m^n: even) 



Pmn = ^ ffm (6*; ) siu 9^1 + 27r g„t (6*; ) < 
1=1 1=1 ^ 

n o /n \ f Ci2 (to: odd,rt: even) 

^—^ [C2i COS t'; (to: even,n: odd) 



(37) 



{-C2i){l3i- P2) (to, n: odd) 
(m,n: even) 



^ [ sm^ 9i sm fc*; J ^ [ C22 m — 

i\T n(R cos9ism9ni A , . / C22(/3i - /32) (m: odd,n: even) 

Af™„ = „2(/3i -/32)^5„(0O— ^^2^ +47r^5™(eO | in i±£2||i (,„: even,n: odd) (^^^ 

In what follows, we adopt 5m(^) = 2 sin [TO0/(iir +1)] as was done so for the 't Hooft model fig] ]. 



C. Mass inequalities 

Since we have at hand the methods for obtaining the physical properties of meson states, we are in a position to 
compute the mass inequalities. For investigating the mass inequalities, we use the properties of the lightest meson 
state in each channel. In Fig. |^ we first plot the behavior of the mass inequality Sfiab/ fJ-ab defined in Eq. (|l|) for finite 
mass differences for a typical case ofG' = G5 = l, (3 = 1. The relative quark mass difference parameter |A| < 1 by 
definition and the mass difference is symmetric with respect to the interchange A —A. At the same time, we also 
plot the behavior expected from the susceptibility, TZA"^. We see that the susceptibility describes the mass difference 
quite well unless the quark mass difference is quite large, say A ^ 0.4. 
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FIG. 3. The normalized meson mass difTerence as a function of A (solid). The mass difference expected from the susceptibility 
(dashes) is also shown. 



Let us move on to the behavior of meson mass susceptibilities. To compute the susceptibihties, we may just use the 
methods explained in the previous section and obtain the susceptibility as the limiting case of small mass differences 
going to zero. While this is logically fine, it incurs unnecessary numerical errors during the process. Therefore, we 
can refine the method by perturbing in the mass difference analytically and then obtaining the mass susceptibilities 
directly. However, the standard perturbation formulas are not applicable to either of the two methods explained in 
the previous section since we are dealing with a perturbation that changes the boundary conditions, as we can see 
from Eq. (jlj). While the formulas that need to be derived should be of use to further study, since this is technical 
and somewhat involved, we have chosen to describe the methods concretely in the appendix. We have computed the 
susceptibility using both methods and have checked that the results do agree. 

The parameters of the gauged four fermi models are /3, G, and the ratios of constituent masses. We expect G 
to play a not so dominant role in determining the properties of the lightest meson state. G5 is the pseudo-scalar 
coupling that strongly affects the lightest meson. f3 is effectively the inverse of the strength of the gauge coupling. 

We first investigate the behavior of TZ with respect to /3 as in Fig. ^. When G5 7^ 0, for large /3, the susceptibilities 
approach those of the generalized Gross-Neveu model, which is quite natural since the gauge coupling is effectively 
weak. This behavior is quite visible for (G, G5) = (0, 1), (1, 1) cases in Fig. ^ and the approach already occurs for 
moderate /3 values, f3 ^ 0.1. When G5 = 0, as we can see from Eq. (^), 7?. = in the generalized Gross-Neveu model. 
In the gauged four fermi model, TZ behaves as ~ f3~'^^^, when G5 = and large /3 as we can see for (G, G5) = (0,0), (1, 0) 
cases in Fig. ^ This is consistent with the expectation from the quantum mechanics calculation in Eq. ( |1C| ) for the 
linear confining potential. For G — G5 and /? = 0, it can be shown that TZ ^ as This behavior is indeed seen 
in Fig. I for (G, G5) = (0, 0), (1, 1) cases. 
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FIG. 4. The behavior of TZ against jS for (0,6*5) = (0,0) (□), (1, 1) (Q), (0, 1) (A), and (1,0) (v)- Dashes indicate /^"^/^ 
behavior and dots indicate behavior. 



Let us now analyze how TZ behaves with respect to G5 as in Fig. |5|. It can be seen that for fixed f3, TZ approaches 
the generalized Gross-Neveu model value as we increase G or G5. Qualitatively, this can be understood as the 
gauge coupling becoming relatively less important when the other couplings are strong. For small G5, the behavior 
is governed by the gauge coupling and we see in Fig. || that the susceptibilities for the same /3 value approach each 
other. While these behaviors can be understood from the physics of the model as we did so here, it is quite non-trivial 
derive them analytically. 
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FIG. 5. The behavior of 7^ against G5 for (G, /3) = (0, 0.1) (□), (0,1) (Q), (1,0.1) (A) and (1, 1) (v)- For comparison, 7^ 
for the generahzed Gross-Neveu model is also displayed for G = (solid) and G = 1 (dashes). 

We have investigated the susceptibility extensively within the parameter space of the theory and found that it is 
positive, except for a relatively small region, which we now discuss. In most regions in the parameter space, the 
numerical convergence of the susceptibility is quite rapid at least in one of the methods and both methods yield 
consistent results. In all these regions, the susceptibility parameter satisfies TZ > 0. However, for small (3, the 
convergence is rather slow. Particularly intriguing is the region G^Gs, 0^/3<Cl. A simple argument shows why 
the behavior in this region can be subtle: In general, the finite dimensional numerical results are analytic with respect 
to the parameters G, G5, /3. From the behavior 7?. ^ as /3 ^ 0, G = G5, we know that unless TZ vanishes as [G—G^Y 
or some higher even power for (3 ~ Q, there will be a region where TZ is negative. Indeed, we find in the numerical 
results that TZ is negative in the region G>G5,0</3^ 1 using both the variational and the Multhopp's method. 
Even the extrapolated values, in some cases, are negative. The meson mass squared is always positive even in these 
cases and the physics of the system seems to be quite consistent. Naively, we would claim that the susceptibility 
and hence also the mass inequality is negative in this regime. This, to our knowledge, would not confiict with any 
general arguments regarding mass inequalities. However, to conclude this would be somewhat premature, since if we 
study the negative region in the parameter space, we find as in Fig. that it shrinks when the basis space is enlarged 
and the negative region is quite small compared to 0(1) which is the natural scale in the problem. It should also 
be noted that even liTZ = ai(G — G5) + 0((G — Gs)^) for /3 = when the basis space is finite, it is still possible 
that the coefficient ai approaches in the full basis space, so that the positive susceptibility is compatible with a 
negative one in the truncated basis space. On the other hand, the regions of negative susceptibility have a common 
region with respect to both methods so it is also possible that a finite region of negative susceptibility remains even 
when the basis space is complete. We therefore conclude that while the susceptibility may be negative in the regime 
G^Gs, 0^/3^1, further investigation is necessary to clarify this point. An analytic computation determining the 
sign of the inequality would be ideal. If this is not possible, a set of basis optimized for the gauged four fermi models 
in this particular parameter regime, in either variational or Multhopp's method should settle this issue. 
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FIG. 6. The zeros of the mass susceptibihty parameter in the G^-jS plane when G = 0.1, computed using finite dimensional 
basis spaces. The solid curves, from top to bottom, correspond to the zeros in the variational method for the basis space 
dimensions of 8, 10, 12, 14, 16, 18, 20. The dashed curves, from top to bottom, correspond to the zeros of the susceptibility 
in the Multhopp's method for the basis space dimensions of 20, 40, 100, 200, 400, 800, 1000. In the small regions below the 
respective curves, the susceptibility is negative. The negative region becomes smaller with the increase in the size of the basis 
space in both methods. 



VI. SUMMARY AND DISCUSSIONS 



We have systematically and quantitatively studied the mass inequalities in gauged four fermi models in (1 + 1)- 
dimensions. Even in the cases where the mass inequality has been shown to be positive from general arguments, 
the size of the inequality is unknown unless an explicit computation is made. We believe that the mass inequality 
is an interesting dynamical quantity characterizing the spectrum of relativistic quantum field theories. To analyze 
the inequalities quantitatively, we adopted a natural susceptibility parameter to compare the size of the inequalities 
throughout the field theory space. We found that the parameter captures the essence of the mass inequality when 
the constituent mass differences are not too large. In the family of generalized Gross-Neveu models, we were able to 
derive an analytic expression for the meson mass susceptibility. In the more general case of gauged four fermi models, 
we have developed methods for obtaining the mass inequalities systematically and have computed them. Since not 
much seems to be known about the quantitative behavior of mass inequalities, we think that it is significant to have 
a class of relativistic field theory models where it has been studied explicitly. While the results are interesting, there 
remain further questions which should be answered. 

An important question is whether the positivity of the mass inequality is much more general than the cases where 
it has been shown to hold j^]. In particular, an intriguing problem is whether there is a relativistic quantum field 
theory wherein the mass inequality is negative yet its physics behavior is consistent. We have found that the meson 
mass susceptibility is positive for most of the parameter space in the gauged four fermi models and have explained 
some of the behavior analytically. The models we studied here include the celebrated models of 't Hooft and of Gross 
and Neveu for special choices of the parameters. For the 't Hooft model, the standard arguments do apply and 
we may show analytically that the inequality is positive. However, in general, no such arguments can be applied to 
gauged four fermi models. Furthermore, the Gross-Neveu models are known to be equivalent to models with Yukawa 
couplings and also display dynamical chiral symmetry breaking behavior. These are exactly the kind of situations in 
which we might doubt that the mass inequality is positive 1^,^. It is interesting that even in these cases, the mass 
inequality is positive, so that in fact, the property holds in much more general than those situations where it has been 
proven. It would be interesting to find an analytic proof for this property if possible. It is important to understand 
why the inequality is positive for the gauged four fermi models and clarify if this can be extended to other theories, 
such as supersymmetric theories with bound states. There still remains a small region within the field theory space 
wherein the sign of the susceptibility, hence also the mass inequality, remains uncertain and further investigation is 
necessary to establish its sign. While (1 + l)-dimensional theories such as the 't Hooft model or the Gross-Neveu 
model have physical behavior resembling those of higher dimensions, we should mention the possibility that in higher 
dimensions, the behavior of the susceptibility might be quite different. Also, even in (1 + 1) -dimensions, the mass 
inequalities might behave qualitatively differently for other class of models. 

We have seen in section |l| that the mass inequality is satisfied in a large class of quantum mechanics models [D . 
This leads us to suspect that the mass inequality is valid for a large class of relativistic quantum field theories also. This 
is certainly consistent with our findings here. However, it should be noted that spontaneous breaking of symmetries 
is essentially a field theoretical behavior, which is also quite relevant to the theory we studied. Therefore, we believe 
that it would be worthwhile to perform further research and in particular, clarify whether the mass inequality can 
become negative in relativistic quantum field theories. In another direction, large N limit of field theories, such as the 
class of models we study, are presumably described by some kind of string theories |^ — an idea, which has recently 
been made more concrete It would be interesting to find out what kind of string theories our models correspond 
to and to elucidate how mass inequalities fit into the string picture. 



APPENDIX A: PERTURBATION THEORY IN THE MASS DIFFERENCES FOR THE SPECTRUM 



Here, we shall briefly outline how to perform perturbation with respect to the relative constituent mass difference, 
A, in the methods explained in section VB for obtaining the spectrum. The standard perturbation methods can not 
be applied here. One major reason is that the boundary conditions ( |l5| ) depend on the masses of the constituents so 
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that they need to be perturbed also. There are additional complications for both the methods used in section VB, 
as we shall describe below. 

In both cases, we perturb in the relative mass difference A and obtain an expansion for the meson mass in terms 
of A, for the cases (Mf, Af|) = {Ml Ml), (M^, M^), (Mf, M^). 

= fil + Afij + AV2 + O(A^) (Al) 

In the first two cases, the first order term exists and are of the same size but of opposite sign and in the last case 
the first order term is absent. Therefore, the leading order term in the mass difference Sfiab will be of order A^, as it 
should be. 



a. Variational method 



In the variational method, we need to consider a generalized eigenvalue problem with the normalization matrix 
not being the identity matrix. In theory, we can just orthonormalize the basis vectors, but in practice, this is not 
numerically equivalent since the normalization matrices can become almost singular even though we have tried to 
normalize the matrix elements to be of order one. Furthermore, since the boundary conditions also are perturbed, 
the normalization matrices will also have a non-trivial expansion in A. 

Let us expand the matrices as 

H = Ho + AHi+A'^H2 + 0{A^), N = No + ANi + A'^N2 + 0{A^) (A2) 
Assume that we have the complete eigen system for the 0-th order problem: 

HoWOn = nl„NoWOn, {wom, NoWQn) = Smn (A3) 

Then, we obtain the expansion for mass squared of the meson state labeled by n 

lAn = (w^On, {H2 - M0n^2) WOn) ~ {wQn, NiWQn) {wOn, {Hi - nl^Ni) WOn) (A4) 
+ 2 ^ 2 K^On, [Hi - Mo«^l) Wom)^ 

We need the expansions of the matrices H,N in terms of A for the three cases, (M^,M|) = 
(M^, M^), (M^, M^), (M^, M^), to obtain the final results. Since this expansion is cumbersome but logically straight- 
forward, it will not be explicitly presented here to save space. 



h. Multhopp's method 

In Multhopp's method, the matrices are not Hermitean so that we need to perform the perturbation theory with 
some care. Furthermore, due to the perturbation in the boundary conditions, the matrix P will also be perturbed. 
To perform the expansion, we will reduce the equation to a mathematically equivalent problem, 

{n"^ - P-^M) w = (A5) 

Ideally, it is better not to invert matrices numerically, but it is a substantially more complicated numerical task to 
solve a generalized non-symmetric eigenvalue problem and also, in this case, the matrix P turns out to be quite robust 
against inversion even for moderately large basis spaces with dimensions of order 10'^. 
We expand the matrices as 

P = Po + APi + A^Ps, M = Mo + AMi + A^Afa 
(P-^M) = (P-Hl)^ + A {P-^I)^ + A^ (P^^Af)^ + 0{A^) (A6) 



where 



{P-'M)^ - Pq-IA/o, {P^'M)^ = Pq-i (Afi - PiPo-'Mo) 

{P-'M)^ = Pq-i (Af2 - P2Po'Mo - PiPo^'Mi + PiPo-'PiPo'Mo) (A7) 
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We need to first solve the 0~th order problem for the left and right eigenvectors, and {wn}Q 

UOm (P^^M)^ = ^l^UQrri, (P^^M)^Von = Mon'^On, (U0m,f0n) = ^mn (A8) 

Then, we may obtain the expansion for the meson mass squared of the meson state labeled by n as 

tJ,j^={uon:iP'^M)^Von) (A9) 
2 r ir. l.r\ \ ^ {uon,{P-^M).Vok) {uok,{P-^M) von) , ^ ^ 

The rest proceeds as in the variational method case. As in the case of the variational method, the explicit expressions 
for the matrices are not shown here due to space considerations. 
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